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Abstract 

1} ■ Supersymmetric t-J Gaudin models with both periodic and open boundary 

conditions are constructed and diagonalized by means of the algebraic Bethe 
ansatz method. Off-shell Bethe ansatz equations of the Gaudin systems are de- 
rived, and used to construct and solve the KZ equations associated with sZ(2|l)W 
superalgebra. 

GO 

I introduction 

In the study of one dimensional long-range interacting systems, Gaudin type models 
ocupied an important place, due to their role in establishing the integrability of 
\q ', the Seiberg-Witten theory [Q, |3| and diagonalizing the BCS hamiltonian of ultrasmall 
metallic grains |4], || |6|. They also served as a testing ground for ideas such as the 
■ functional Bethe ansatz and general procedure of separation of variables [J7|, [5], [|. 

The t-J model was proposed in an attempt to understand high-T c superconductivity 
TfJ| . It is a correlated electron system with nearest-neighbor hopping (t) and anti 



O 
O 

(N 

X> 

(N 
(N 



ferromagnetic exchange (J) of electrons. Using the nested algebraic Bethe ansatz 



method, Essler and Korepin obtained the eigenvalues of the periodic system |TT| . Soon 
after the open boundary case was studied in [13, 113 . 



In [14], the periodic t-J Gaudin model was investigated and its eigenvalues were 



^ . obtained. In this paper, we study both the periodic and open boundary t-J Gaudin 
models by a method different from that used in ||14j| . 

The Knizhnik-Zamolodchikov (KZ) equations were first proposed as a set of differ- 
ential equations satisfied by correlation functions of the Wess-Zumino-Witten models 
||15|| . The connection between Gaudin type magnets and the KZ equations has been 
studied by many authors ]Tj| [H], [TS|, [19], [20| . We are interested in the super KZ equa- 
tions associated with sZ(2|l)W superalgebra. We will construct and solve these KZ 
equations with the help of the t-J Gaudin models. 

The outline of this paper is as follows: In the first part of this paper (section 2- 
4), we study the periodic t-J Gaudin model and its corresponding KZ equation. The 
t-J model hamiltonian is constructed in section 2, and diagonalized in section 3, by 
using the algebraic Bethe ansatz method. We also derive its off-shelled Bethe ansatz 
equations, and use them, in section 4, to construct solution to the corresponding KZ 
equation. The second part of this paper is devoted to dealing with the open boundary 
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t-J Gaudin model. In sections 5 and 6, we construct and diagonalize the open boundary 
t-J Gaudin model. In section 7, we obtain solution to the KZ equation associated with 
the boundary Gaudin model. 



II Periodic t-J Gaudin model 

II. 1 t-J model with periodic boundary condition 

The supersymmetric t-J model is described by the R-matrix arising from the 3-dimensional 
representation of sZ(2|l) < - 1 ' ) . In the Fermionic, Fermionic and Bosonic (FFB) grading of 
the representation space, the R-matrix is given by plj 
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(11.1) 



where rj is a crossing parameter and 



a(A) = 
c + (A) 



6(A) 
e A sinh(^) 



sinh(A) 
sinh(A — if 



w(X) 



sinh(A + rj) 



c_(A) 



sinh(A — t]) ' 
i~ x sinh(^) 



sinh(A — 77) ' sinh(A — if 

This i?-matrix satisfies the graded Yang-Baxter equation (YBE) 



(II.2) 



(11.3) 

where e a is the Grassman parity: e a = for bosons and e a = 1 for fermions. The 
R-matrix satisfies the unitarity and cross-unitarity relations, 

i?i 2 (A)-R 2 i( — A) = p(A) • id, p(A) = — sinh(A + rf) sinh(A — rf), 

Rfi{X - r?)Mii?2i ri Mf 1 = p(A) ■ id, p(X) = sinh(A) sinh(A - rf), (II.4) 

where M is a diagonal matrix diag(e 2,? , 1, 1) and st is the super-transposition defined 
by 

(A% 3 = A^-lp^ . (II.5) 

Consider the L-operator 



L aq (X) = R aq (X), 



(II.6) 
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where a represents the auxiliary space and q represents the quantum space. The L- 
operator also obeys the (graded) YBE 

R 12 (X - /j,)L 1 (\)L 2 (h) = L 2 ( f i)L 1 (X)R 12 (X - (x). (II.7) 
The tensor product is graded, namely, 

{F®G) b a d c = F b a G d c {-)^ + ^. (II.8) 
The row-to-row monodromy matrix Tjv(A) is defined as the product of N operators, 

T a (A) = L al (X - Zl )L a2 {\ -z 2 )--- L aN {X - z N ), (11.9) 

In matrix form, 



1"JV 



(11.10) 

By repeatedly using the YBE, one can easily check that the monodromy matrix satisfies 

R(X - //)Ti(A)T 2 (/i) = T 2 (/i)T 1 (A) J R(A - /i). (11.11) 

The transfer matrix t(X) is defined as the supertrace of the monodromy matrix over 
the auxiliary space: 

t(A) = strT(\) = ^(-l) e «T(A) aa . (11.12) 

Using the YBE, one can show that the transfer matrix t(X) constitutes a one-parameter 
commuting family, i.e. 

[t(\),t(ji)] = 0. (11.13) 
Therefore, the t-J model is integrable. 

II. 2 Supersymmetric t-J Gaudin model 

Supersymmetric t-J Gaudin model can be obtained by taking the quasi-classical limit 
7] — > of the transfer matrix at the point A = Zj |[22|| . So we expand the R- matrix, 
L-operator and transfer matrix around the point rj = to get 

R(X) = l+7 1 R(X)+0(r ] 2 ), (11.14) 
L(A) = 1 + rjC(X) + 0(r] 2 ), (11.15) 
t( Zj ) = -1 + r)i{ Zj ) + Oirf). (11.16) 

Then the Hamiltonian of periodic t-J Gaudin model can be obtained from the second 
term of (PXToT). 
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From (|II.14|) , we have 

Rij(X) = 2 coth(A)e 33 <g> e 33 + coth(A) 



+e 



22 



33 



e ll + e 33 



e n ® e 22 + e u d$ e 



e 22 



3 

33 



+ e" 



sinh(A) 
1 

sinh(A) 



e 13 ® e 31 + e 23 ® e 32 e 12 <8> e 21 



'31 



e 13 + e 32 ® e 23 — e 21 ® C 12 



(11.17) 



5j n 5jfi. Denote 



where ejy is a matrix acting on the k-th space with elements (e^)^ = S ia Sj/3 
by S 1 , S^, S" 2 , <5±i, Q±i and T 2 the generators of s/(2|l), which satisfy, among others, 

[S\S] = S Z , {Q\,Q 1 } = T Z , {QU,Q-i} = S z + T z . (11.18) 

In the fundamental representation, they take the form, 

I 1 °\ 

S z = 0—10, T z = 

V o o o / 

S = e 2 i, 5 f = ei 2 , Qi = e 32 , 

<2i = e 23 , g_i = e 31 , Ql^eia. (IL19) 

Thus, 

4(A) = 2coth(A)[(l-(^) 2 )(l-(^) 2 ); 

+ coth(A) [(1 - T z ){2 -T z -S z - (S z ) 2 ) + (2 - 2? - S 2 - (S 2 ) 2 )(l - 2J) 




+ (1 - 2? - 5?)(1 - + (1 - (S z Y)(l ~T Z -S 



+ 



sinh(A) 



-SjSj + ^ Qi jC rQi,a 

<T=U 



+ - 



sinh(A) 



SiSj + ^ Qi,aQi,a 



(11.20) 



Using the standard fermionic representation [11 



Sj — Cj^c^—i, Sj — Cj_iCj } i, Sj — n j,i> 



rjnz 



1 - 7l 3 -,_i, Q} i(T = ct )(J (l - 7lj- _ CT ), Qj !<T = c j!lT (l - %- CT ), (11.21) 



where n^± = (i ± Cj t ± and rij i+ + n,j_ = rij, we have 

4(A) = 2coth(A) [(1 - (n u - ^, T ) 2 )(1 - (n i4 - n i)T ) 2 )' 

+ coth(A) [n u (l + n j>T - (n jd - n iiT ) 2 ) + (1 + n iiT - (n^ - n ijT ) 2 )n ia 
+fii, T (l - (fii,j - n iiT ) 2 ) + (1 - (n i4 - n iA -) 2 )n jA 
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+ - 



sinh(A) 

A 



+ — 



-4,T c ia c Ja c i>T + E Ci,«r(l - rai,- ff )ct ff (l - 
<r=T4 



sinh(A) 

Similarly, from the eq.( |11.15| ), we obtain the elements of the Lax matrix 

£i X = coth(A) — coth(A)nj _i, C 3 2 2 = coth(A) — coth(A)n Ji i, 

£33 = coth(A) + coth(A) 1 — (iij-i — %i) 2 , 
.-a p \ 

+ rJ + 

smh(A) • y ' smh(A) J ' 

„— A „A 



(11.22) 



13 



sinh(A) 



^23 



smh(A) J ' 

Ci.iC 1 - £32 = ^7TT c A-i( 1 ~ n i.0> 



H 



sinh(A) ^ sinh(A) 3 

Finally, the Hamiltonian of the t-J Gaudin model can be written as 
dt(X = Zj) 



(11.23) 



drj 



N 



|?7=0- 
1 



V - 

k=x+j smh(z j - z k ) 



+ e 



-(Zj—Zk) 



{cosh(z i - Zk)[—Tij-ink,-i - n jt ink,i + (1 - nj)(l - n k ) 
c] a (l - nj- a )c k ^{l - n k - a ) - S]S k 



<T=±1 



+ e 



Z j—Zk 



c,>(l - n i ,_ <r )4 ))T (l - n fc _ CT ) - S^J 

T=±l 

III Bethe ansatz for the t-J Gaudin model 



(11.24) 



Hamiltonian ( 11.24 ) can be diagonalized by using the algebraic Bethe ansatz method. 

(III.l) 



To simplify calculation, we use the gauge transformation 
R 

to get 



R{\) 
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(e A / 2 ,e A / 2 
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where d(X) = sinh(?y)/ sinh(A — rf). The corresponding L-matrix can be written as 
L n (X) 

6(A) - (6(A) - a(A)K x -C-(A)e^ d(X)e^ 

-c + (A)e™ 2 6(A) - (6(A) - a(A))e™ 2 d(X)e™ 2 

d(X)e n 13 d(X)e n 23 " 6(A)- (6(A) - w(\))e% 3 

(III.3) 

We expand the row-to-row monodromy matrix (|II.9|) around r\ = 0: 



T(A) 



f A 11 (\) A 12 (X) B 1 {X) 
A 21 (X) A 22 (X) B 2 (X) 
\ Ct(X) C 2 (A) -D(A) 



1 + rjf (A) + e>(7? 2 ) 
l + ?7 



f i n (A) Ai a (A) 4(A) \ 

i 2 i(A) i 22 (A) 5 2 (A) 

\ C 1 {X) C 2 (X) D(X) J 

By the graded YBE (|TTTTT1) , one finds 

[fi(A),f 2 0i)] = [Ti(A) + T 2 (fi), R 12 (X - //)] 

and the commutation relations, 

C Sl (ni)C S2 (^ 2 ) = -C S2 (n 2 )C sl (nA, 



(III.4) 

(III.5) 
(III.6) 



D{zj)C s (n) 



coth(z,- - n)C s (/j,)D(z j )\, n=0 + C s (^)D{zj 



sinh(/i — Zj) 



sinh(/i — Zj 



-C s (zj)D{n) v= Q, 



(III.7) 



+ 



r MN( Z j ~ / i )) r? =0 C*s(/ i ) y 4sim(%)-'?=0 + C pi (fx) A SlS2 (Zj 

~- w ~ \C S2 ( z j)A SlPl (/i) r?= o, 



(111.8) 



where j indicates the lattice position and Ej (s) acts on the quantum space with 
Ej(s) = ei 3 for s = 1 and E~(s) = e 23 for s = 2; the r-matrix r(A) is defined by 



r(A) 



sinh(A) 



^ sinh(A — 77) \ 

sinh(A) — e~ A sinh(7/) 

— e A sinh(?7) sinh(A) 

V sinh(A-77) / 



(III.9) 
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Define the vacuum state: 



10 >n- 



( 



V i 



10 >= 



i|0 > fc . 



(111.10) 



Then, 



5 a (A)|0>=0, C a (A)|0>^0, 



N 



i=l 



where j = 1,2,- 

Define the Bethe state 



D(A)|0 >= ^2coth(A - zi)\0 > 

A ab (X)\0> 
■ ■ ■ , N. 



for A = Zj, and a ^ b 
Y^iLi coth(A — Zi) for \ Zj, and a = b 



(in.ii) 



(111.12) 



= CaMC^te) ■ --CdMlO > F dl "' dn , 
where F dl " dn is a function of the spectral parameters fij. We moreover define state 

^ = C^\^)C^\^)...C^\^). (111.13) 

Then 0« spans a subspace of the space spanned by 0. Applying the transfer matrix 
to these states and keeping in mind of the commutation relation ( |111.7| ), we can find 
the eigenvalues Ej of t(zj) and the Bethe ansatz equations. This is done as follows. 
Firstly, from ( [[11.11 ) we have 



N 



coth(zj — Ha) — 2 coth(zj — Zk) 



a=l 
n 



i) 



a-l 



\ sinh(/i a — Zj) 



N 



2 coth(// a - z k ) 

. k=l 



E, (d a )(f) c 



- coih(/i a - up) 

n l 1) Q_1 



where Q = nja=i,=t a Cdpiflp)- Secondly, the action of A(zj) on <p is given by 



(111.14) 



A-aai^Zj , 



E- 



\a-l 



j sinh(/i a — z 



N 



+ E 



sinh(/x a - Zj 



Y COth(/i Q - Zfc) + A (1) (/i a ) £y (d a )0 Q 
J-' L fe=l 
1)«-1 . 

Cd a {Zj)(f>a, 



(111.15) 
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where AW(^ a ) is the operator 
d 



f«(A) 



sir 



r?=0 



Using results in the appendix, we obtain 



A (1) (/ia) = J] coth(/i a - £t/?) - E coth (^c - 

P=l,^a 7=1 



(111.16) 



(111.17) 



where /4 , ■ • • , satisfy the constraints 



/«= £ coth(^ 2 £ coth(/iW - = 0. 



5=1,^7 



(111.18) 



Thus, from (|L2g) ,( |lLTg) -( glLTg ) and ( |ILTg )- |lLTg ), we obtain the off-shell 
Bethe ansatz equations 



M 



E ^ + E — T77 1 —f<*Ej{d, 



a ^ sinh(/i Q - zj 
where \x a satisfy the condition = and 



(111.19) 



E; 



fa 



ri N 

coth(zj — fj, a ) — 2 coth(2; J - — z k ) 

a=l k=l,^j 
m N 

- J2 coth(/x a - fijp) + coth(/i a - z k ), 

7=1 k=l 
n ( N 1 \ 

II E — T7 lo > i^-** 

j-Ji V^i smh(/i a - Zfc ) ; 

n / JV i \ 

n E-T7 — ^{d p ) iof 4 "* 



IV Super KZ equation 

As a set of partial differential equations, the KZ equations take the form 

= for j = 1,2,- --,N, 
where the differential operator Vj is defined by the Gaudin Hamiltonian Hj\ 



(111.20) 
(111.21) 
(111.22) 
(111.23) 



(IV. 1) 



(IV.2) 
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with k being a parameter. Substituting (|II.24|) into ( IV.2|) , we can check 

[V„V fc ]=0, (IV.3) 

which ensures the integrability of the KZ equation. 

To simplify our calculation, we make the following transformation 

JV 

Hj — > Hj + 2 coth(zj — Zk), 

N n 

Ej — > Ej + 2 coth(zj — 2fc) = 2J coth(zj — fi a ). 

k=l,^j a=l 

Under this transformation, the form of the off-shell Bethe ansatz equations is invariant. 

The function *&(z) can be constructed by the hypergeometric function x(z, fx) which 
oebys the equations 

d 



K dz~ X = E ^ 



d 



faXi 



and the constraint = 0. The solution to the above equations is given by 



n N 

X(z,») = n [sinh(// Q - ^T l/K II Il[ sinh (^ - t*a)] l/K 

(3<a a=l j=l 



(IV.4) 



(IV.5) 



with \i a satisfying the condition = 0. With the help of x(z,fi), the function ^(z) 
is given by 

=f"-f d »i ■■ ■ df*nX(t, z)<Kt, z), (IV.6) 

where the integration path C is a closed contour in the Riemann surface such that the 
integrand resumes its initial value after t a has described it. Substituting the expressions 
of Vj and *&(z) into (|IV.1| ), we can show that the KZ equation is satisfied. The proof 
is as follows 



9 ,, V 



/ • ■ • I dfii---dfin \K^(j> + KX^~ 
JC JC \ OZj OZj 

j> c ■ ■ ■ j> c dh ■ ■ ■ d/i n (^xEj4> + K *^j 



f ■ ■ ■ f dfxi ■ ■ ■ da r 

JC JC 



xHj(p - xE W(^«> z j )f a E j (d a ) 



a=l 







- — [W(^ ol) Zj)(j) a E j (d a ) 

j> ■ ■ - <j>^ dux • ■ • d[i n [xHj(p 



d 

— (xW(Ha, Zj)(j) a Ej(d c 



d/jL. 



(IV.7) 
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where W(fx a ,z) = (—1)" / sinh(/z a — z). 

V Open boundary t-J Gaudin model 

In this and the next section, we discuss the open boundary t-J Gaudin system. We 
start with the graded reflection relation fl"2|| , 

A(A - + ^)gf 1 1 ^(/^)S(-l) (ebl+eCl)eb2 

= KM%R(\ + f,) b a \llK(\)tlR(\ - AiJSJC-l)^ 4 ^^, (V.l) 

where the K(X) is the reflection K-matrix. The diagonal solutions of the reflection 
equation were found in [I2|. In the present paper, we only consider a special case in 
which K(X) = 1. 

Following the standard procedure, we define the double-row monodromy matrix 

T(A) = T{X)K(X)T- 1 {-\). (V.2) 

Here T(A) is same as in the periodic case. One can check that the following relation is 
satisfied: 

= TQi)%R{\ + ii)i&T{\)%R{\ - //)f^(-l)^ +£ -^. (V.3) 

The dual reflection relation reads 

R X2 {H - \)K+(\)M{ l R 21 (r] - A - ^)K+{^)M 2 l 

= K+{n)M^R n (r] - A - fJl )K+(\)Mr 1 R 2 i{v - A), (V.4) 

where M is a diagonal matrix M = diag(e 2r? , 1, 1). Solution K + associated with K{\) = 
1 is given by 

K + (\) = MK{-\ + r}/2) = diag(e 2,? , 1, 1). (V.5) 
Define the boundary transfer matrix, 

t b {\) = strK + (\)T(\). (V.6) 

By ( [11. 4| ) and ( |V.3| ), one can show the commutativity of the transfer matrix for different 
A values. 

Similar to the periodic case, the boundary t-J Gaudin system can be obtained by 
expanding the boundary transfer matrix at the point A = Zj around rj = 0: 

t b (X = zj) = 1 + vHj + Oirf). (V.7) 

The second term on the right hand side gives the Hamiltonian of the open boundary 
t-J Gaudin model. Explicitly, 

_ dtjzj) 
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2(1 + 3coth(2zj)K 



A/ 



+ J2 a ^u, _ „ ^ { C0S H z j - Zk)[-nj,-Wk,-i ~ n jA n k>1 + (1 - nj )(l - n k ) 

k=l=£j bnm \ z j z k) 



+ e 



Zj Z\~ 



■=±i 

C oA l ~ %-^)4,a( 1 - n k,—o) ~ SjSl 



CT=±1 



N 



+ E - 



fc=1 ^. sinh^- + z k ) 



{cosh(z 3 - + 2 fc )[[-n j _in fc _i - n jtl n k)1 + (1 - n^l - n*.)] 



-=±1 

^ c] i0 .(l - rij _ CT )c fej(J (l - rik-o) + S]Sj 



(V.f 



VI Bethe ansatz for boundary t-J Gaudin model 



As in the periodic case, we write the double-monodromy matrix as 

T(A) 



( A n (\) A l2 {\) B 1 (X) 
A 21 (X) A 22 (X) B 2 (X) 
\ C X (A) C 2 (A) P(A) 



(VI. 1) 



Around 7/ = 0: 



T(A) = 1 + r/T(A) + 0(rf) 

( in (A) ii 2 (A) Bi(A) \ 



1 + 77 



^21 (A) ^2 (A) Ba(A) 
V C X (A) C 2 (A) P(A) j 



o(v 2 



(VI.2) 



Applying T(A = Zj), j = 1, 2, • • • , N, to the vacuum state ( |ill.l(J| ), we have 
4(%)|0>=0, (7 a (^)|0>^0 

iV 

£>(z 3 -)|0 >= ^ 2(coth(z j - Zi) + coth(^- - Zi))\0 > . 



i=i 



4 mm — / ^ for A = 2,-, and a ^ b 

A ab (A)\U >= < j2? =1 (coth(\ + Zi ) + coth(A - Zi)) for A ^ z h and a = b ' 



The Bethe state of the boundary system can still be taken as 

b = c d Mc d M ■ ■ ■ c dn (fi n )\o > F d ^- d ". 



(VI.3) 
(VI.4) 
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Write 



i(AU|r,=0 = ^(AU|„=o + $ab . 7^rrr V(\) v=0 , (VI.5) 

smh(2A) 



where 

N r 

A^U\0 >= £(coth(A + + coth(A - Zi )) - — ^— V{X) V=0 . (VI.6) 
i=1 smri^zAj 



Then 



= -A aa {zj) +V{zj) 



K) v=Q A{z 3 ) ri=0 -UV{X) v=Q , (VI.7) 

where a = 1,2 and U = 2/ sinh(2zj). The last term in (|VT.7|) corresponds to the 
boundary condition. 

We now find commutation relations between A a b(\), T>(X) and C^(/i). After a 
tedious but direct computation, we get 

C dl (nx)C d2 (n 2 ) = -C C2 (n 2 )C cl (ni) 7 (VI.8) 

+ . i (--B/(<i)*(M) + C.i(-j)f MI„-o) 

~ ■ u/ 1 i — V (-Ej{d)A bd ([i) +C b (Zj)Abd(n)\r,=0 

SVO.ll[Zj + fl) \ / 

2coth(2u) „ . . . , 2 coshfz, • + w) ^ ,, s 7 / m 

+ ss^^^wi- - sM^+r) £7(6)AiMu ' 

(VI.9) 

Aidi(^')Cd 2 (^) = C d2 (/i)A ldl (^-) 

+ (ria(^ + // + r?)^ 2 r 2 i(^ - A*)**)^ 4(^)^6, (A) |, =0 

+ gin^^. _ ^^ufcAi* ( — -^j (^i)Anda(^) + Cd 1 (z j )A ai d 2 (fi)\n=o 

- ■ w , / aida^Mi ("^(M^G") + C b2 (Zj)V(^)\ n=0 

/sinh(r ? )r 12 (2^ + r ? )S i V / 



sinh^ — //) 



r?=0 



(6l)A 2 d 2 (^)|r?=0 



^ sin(2 / i)^(r / )r 12 (2z J + 7 ? )g 2 ^ 
+ ( sinh(, J + / , + ,)sinh( 2/ , + ,) J^ ^ (VU0) 
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V(zj) v=0 C d (fi) = Qi^V^Zj)^ + . . 1 — r E d (d)V(fi) v=0 

E-{b)A M {v) v =o, (VI. 11) 



sinh^- + /x) 



Aiidi (Zj) v=0 Cd 2 (/i) — Crf 2 (/ij^Laxdx (% )?7=0 + ■ w ~ T-^j (^l)^aid 2 (/ i )??=0 

' S^Ejid^V^)^ (VI. 12) 



sinh(,Zj + /i) 



Then, applying (|VL7|) to the Bethe state and using the above commutation relations 
repeatedly, we obtain the off-shelled Bethe ansatz equations 



H b,b = dt(\ = Zj) 



r)=0S 



dr) 

n 

E)<? - ^ zM'Erfa (VI.13) 



a=l 



where 



^ /£ 2sinh (/ ,)cosh(z fc ) \ 
a=i Vfc=i smh(/i Q - z fc ) smh(/i a + z k ) ') 

€ = ft (E • J sinh ^ )co y fc) ^ - W > )|o>^. (vi.15) 

2 N 
smh(^) fc= f-^. 

n 

+ ^ [coth(z j + /i a ) + coth(zj - /i a )] (VI. 16) 

2 n 

^ = — r Uo — V~ £ (coth(/x Q - ^) + coth(// Q + ^)) 

Smtl^/i tt J /3=l,^a 
N 

+ ^ (coth(/i a - z k ) + coth(// a + z fc )) + A i X) W, (VI.17) 
fc=i 

W 6 (//a,*) = ( ~ ir " l2Sinh( ^ )COSh(/ia) . (VI.18) 

a ' 3 sinh(zj + fi a ) sinh(zj — fi a ) 

Here is the eigenvalue of the nested transfer matrix 

4 1} (A) 

= A [strK«+r(A + A*i + ^ c e V(A + ^ + v )%% ■ ■ ■ r(A + ^ + i?)*^ 

^ (1) r 21 (A - M „)S:- n ld " ■ • • r 21 (A - A* 2 )fr 21 (A - f, n ) b b f e A , (VI.19) 
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where 

K {1)+ = diag (e 2r? , l) , K {1) = (1 - sinh(^) / sinh(2A + 77)) • id. 

The detailed calculation of the nested eigenvalue is discussed in the appendix. Here 
we write the result 

K GO = ^-tt — 7+ [coth(/i a -/^) + coth(/i a + /^)) 

Sinn^/iaJ p=\,^ a 
m 

- ^ [coth(/x Q - + coth(/i a + , (VI.20) 

7=1 

where /A 1 ) satisfies the nested Bethe ansatz equation 

n 

ff = EHM^+^ + coth^ 1 )-^)) 

/3=1 

m 

-2 £ (coth^+^^ + coth^-^)) 

<S=1#7 

= 0. (VI.21) 

Substituting the nested eigenvalue into (|VI.17|) , we obtain the Bethe ansatz equations 
for the boundary t-J Gaudin model 

g-Ma N 

fa = T7 — r + Yl (coth(/x a - z k ) + coth(/i a + z k )) 

cosn(^ a j 

m 

- Y (coth(/i a - + coth(// Q + //«)) 

7=1 

= 0. (VI.22) 

VII Super KZ equation in the boundary case 

As in the periodic case, the KZ equations are 

V i * = 0, V^k-^-H*, j=l,2,..-,N, (VII.l) 

but now Hj is the hamiltonian of the boundary t-J Gaudin model. We make the 
transformation 

N 

Hj — > Hj — 2/ sinh(22 J -) +2 V] [coth(2j + z^) + coth(zj — z k )} , 

N 

E) -> - 2/ sinh(2^) + 2 ^ [cothfo- + z fc ) + cothfo - z k )\ 



[coth(zj + yU a ) + coth^j — /^ a )] 



a=l 
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This transformation leaves invariant the form of the off-shell Bethe ansatz equations. 

To construct ^(z), we introduce a hypergeometric function x( z ^) which satisfies 
the following equations 



K 



d 
d/i a 



(VII.2) 



and the contraint f^ m = 0. Solving these two equations, one gets 



n(i 



a=l 



1/k 



JJ [sinh(/i a - /jty) sinh(/i Q - /i^)] 



<*=1 



(VII.3) 



n AT 

x II Il[ sinh (- 2 i + /■*<*) sinh(^ - Ai a )] 1/K , 

a=l j=l 

where //W satisfies the nested Bethe ansatz equation f^ w = 0. With the help of x(z, /■*), 
the function ^(z) is given by 



^z)= f ■-■]> d f M 1 ---d f M M x(t,z)<f>(t,z), (VII.4) 

where the integration path C is a closed contour in the Riemann surface such that the 
integrand resumes its initial value after t a has described it. Substituting the expressions 
of Vj and *&(z) into fllV.lJ ), we can show that the KZ equation is satisfied. The proof 
is as follows 



8 „ 



■ ■ ■ d/ii ■ ■ ■ dfi M ^ + K * J^~) 
fc " ' fc dtl " ' d ^ M + KX Jz~) 



f ■ ■ ■ f dfii ■ ■ ■ dfi M 
Jc Jc 



M 



xH^ + xY, W (^Zj)f a Er(d 



3 r«yr a 



a=l 







xW^z^Ejida) 



C JC 



d 

, 1 9fi a 



(VII.5) 



where 



W(fi a , Z) 



(-l)«- 1 2cosh(z j )sinh(^ a ) 
sinh(/i a + Zj) sinh( / u a — Zj] 
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A Details on the nested Bethe ansatz 
A.l Periodic case 

The nested Bethe ansatz was used to obtain eigenvalues of the transfer matrix con- 
structed from the following r-matrix: 



r(A) = 



sinh(A) 



/ sinh(A - 77) \ 

sinh(A) -e~ A sinh(r]) 

-e A sinh(r]) sinh(A) 

V sinh(A-r?) J 



(A.6) 

One can check that this r-matrix satisfies the unitarity and cross-unitarity relations 

r 12 (A)r 2 i(— A) = — sinh(A + 77) sinh(A + 77) • id., 

r$(2ri - A)r^(A) = sinh(A) sinh(2A - 77) • id., (A. 7) 

and the YBE 

r 12 (A - fj,)L^\\)L^\fj) = 4 1) ( /J )LS 1) (A)r 12 (A - »), (A.8) 

where 

L W (X ) = 1 ( Sinh ( A " ~ a i Sinh(r?) \ ( A 9) 

3 sinh(A) ^ S inh(77) sinh (A - 77^) J 

with a ± ,a z being the usual Pauli matrices. The monodromy matrix is 
T {1 \\) = Lg\\-^)...lg\\-to)LM(\-vL 1 ) 

= ^ C«(A) £)W(A) J' (A - 10) 

which satisfies 

r 12 (A - ^)T}\X)T}\^) = T 2 1) (/i)T 1 1) (A)r 12 (A - /i). (A. 11) 

The transfer matrix of the nested system is defined by 

t«(A) = strT (1) (A). (A.12) 

One can easily prove the commutativity of the transfer matrix: [^(A), t^(fj)] = 0. 
Expanding r,L^,T^ and T^> around 77 = 0: 

r(A) = l + r 7 f(A)| r?=0 + C9(77 2 ), 

L«(A) = l + r 7 ^ 1 )(A)| r?=0 + C9(77 2 ), 

T«(A) = l + 77f( 1 )(A)| r , =0 + O(77 2 ), 

f«(A) = -2 + r ] P\\)\ ri=0 + O(r ] 2 ), (A.13) 
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where 



f(A) 



/ -coth(A) 









-e A /sinh(A) 






-e~ A /sinh(A) 









-coth(A) J 



4 1} (A) 



T«(A) 



*«(A) 



1+3 

2 

.+„A 



coth(A) -a fc e _A / sinh(A) 



V -o-fc e A /sinh(A) 

/ A«(A) BW(A) \ 
1, CW(A) Z) (1) (A) J 



l-o-f 



coth(A) 



(A.14) 



str \ c-wJa) D«(A) ) 

By the nested YBE ( |A.11| ), we find the following commutation relations 
C^(X)C^(h {1) ) = C^(fi^)C^(X), (A.15) 



i (1) (A)(7 (1) (/U (1) ) 



- coth(A - ^)C^(^) + C'«(//«)iW(A) 

+ sin ^ A A+ _ M ^ (1)) (c (1) (A)i (1) (/^ (1) ) + C«(A)D«(^)) 



coth(A - ^)C^(^\ =0 + C^(^)D^(X) 
+ , — (cW{\)DM(iiM) + C7 (1) (A)i (1) (/i {1) ) 



T)=0 

(A.16) 



sinh(A- / u( 1 )) 



For the nest transfer matrix, we choose its vacuum state as 



\o>[ 



»?=0 

(A.17) 



(A.18) 



Applying the T-operator to the vacuum state, we have 

i (1) (A)|0 > (1) = £«(A)|0 > (1) = 0, C7 (1) (A)|0 > (1) ^ 0, 

m 

D«(A)|0>«=^coth(A-/i,) \^fi k (A; = 1,2, ••• 



m) 



(A.19) 



The eigenvector of the nested transfer matrix ([A. 14 ) can be constructed from 

^ = C^\^)C^\^)...C^ (A.20) 

Applying the nested transfer matrix to the above state and keeping in mind A = fi a , 
we obtain 



A (1) (/i Q ) = coth(/x a - y.p) - coth(/i Q - /jW), 

/3=l,^a 7=1 



(A.21) 
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where 4 , ■ • • , ^ satisfy the constraints 

n m 

/«= £ coth(^-/i«) + 2 £ coth^-^HO. (A.22) 
A. 2 Open boundary case 

In the open boundary case, we have boundary reflection and dual relection equations 

r 12 (A - /x)Kj 1) (A)r 21 (A + ^K^fa) = (fi)r l2 (\ + /i)^ 1) (A)r 21 (A - //), 

(A.23) 

r 12 (/x - A)K 1 (1)+ (A)Mf 1 r 21 (2r ? - A - (ji)M^ 
= (^M^r^r, - A - //)irP + (A)Mf 1 r 2 i(/i - A). (A.24) 

One can check = 2 cosh(?y + A) sinh(A)/ sinh(2A + rj) and i^ 1 )" 1 " = diag(e 2,? , 1) are 
solutions to the first and second equations, respectively. 

Using the nested monodromy matrix ( |A.10| ), we define the double-row monodromy 
matrix for the open boundary system 



T (1) (A) = T (1) (A)ir (1) (A)T (1)_1 (-A) 
I C«(A) DW(A) J ' 



(A.25) 



where T^ 1 ) and 1 are defined by 

T^ n (X)t^Z = r(A + /i 1 )^r(A + /i 2 )2g...r(A + /i B )^ ieii 

= 4 1) (A + /ii)4 1) (A + /i2)---4 1) (A + /in), (A.26) 

= LW _1 (-A + ^)...4 irl (-A + /22)4 irl (-A + /*i), (A.27) 
respectively, and the L-operator takes the form 

r m m - ( & (A) - (6(A) - a(A))e? -C-(A)e| 1 \ 

fe (A) " I -c + (A)e^ 2 6(A) - (6(A) - a(A))ef J " (A ' 28) 



Let A = A + 77/2, jl = jj, — r]/2, one sees that the above definitions coincide with ( |V1.19|) . 
The double-row monodromy matrix satisfies the reflection equation 

r 12 (A - /U )T 1 (1) (A)r 21 (A + ^)T 2 (1) (//) = T 2 (1) (/i)r 12 (A + /U )T 1 (1) (A)r 21 (A - //), (A.29) 

Thus, we can define the transfer matrix as 

4 1} (A) = strif (1)+ T (1) (A). (A.30) 
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Around r] = 0, we have the expansions 

tf>(A) = -2 + 4 1) (A) r , =0 + O(r ? 2 ). (A.32) 

Write 

- (i) p- 2A 

^)(A)| f;=0 = ^ ( X )U=o-^2\f {1) W\v=o- (A.33) 

The we find the commutation relations between t^\\) and (fi^) 

C«(A)C«(/i«) = C«(^)C«(A). (A.34) 

+ ^r^ ^) + C«(A)^))) (A.35) 



^ (1) (A)C (1) (/^ (1) ) 



Sinh(2A) -C^\^)A^\X) v= o + C^\^)A 1 \\) 



sinh(A + sinh(A - 

„-(A+m (1) ) 

._^^( C a»(A)^) W + d(.»(A) I >(-» M )^ 

^ <A_ " < " ) fc< 1 >(A)^ (1) ( M < 1) )+C< 1 >(A)i< 1 >( M < 1 >)) , (A.36) 

/ ?7=0 



h(A- / u( 1 )) V 



sin 



P (1, (A)C(/. (1) ),=o = C(^ (1) )P (1) (A), = o + g . nh e (A _ ^ ^(A^W^Vo 

-C^(X)A^(^) V=0 , (A.37) 



sinh(A - 

^(A)^)^ = C(//)^ 1 )(A) + g J h( ( ^ ^ CM(\)AU(») v= o 

p -(A+/i) 

- ss^) c,,, ( A)I " ,,(ri -»- (A - 38 > 
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Substituting (A) = 2 cosh(r/ + A) sinh(A)/ sinh(2A + rj) and (A) = diag(e 2r? , 1) 
into the nested transfer matrix, one obtains the transfer matrix (A) of the t-J Gaudin 
model 



f(D ( ~ x) = d str (e 2v 0W^)(A) BW(A) 
* 6 W " d v Str { ill CW(A) 2>«(A) 



r)=0 



= _i(i)( A )_p(i)(A)_2^W(A) J7=0 

- (1) p- 2 ^ 

= -.4 (A) - X> (1) (A) - 2A ( - 1 \\) V=0 ——~-V^\X) v=0 . (A.39) 

smh(2A) 

As in the periodic case, define the vacuum state for the nested open boundary 
system 

|0 > (1) = <g>JJ =1 |0 >1 1} (A.40) 
with |0 >j^ = ^ ^ ^ • Then, applying the elements of T(A) to the vacuum, we obtain 

B«(A)|0>=0, C (1) (A)|0 >^ 0, (A.41) 

P«(A)|0> = -( —!— + f>oth(A + ^) + coth(A - ) |0 > (1) 
\sinh(2A) ^ / 

A^ fc (fc = l ) 2,-.,m). (A.42) 



Using the nested YBE (|A.11|) and the transformation ( |A.33| ), we obtain 



1 (1) (A)|0 > (1) = - A, ttt |0 > (1) • (A.43) 
smh(A) 

The eigenvalues of the nested open boundary Gaudin system can be obtained by 
applying (X) to the eigenvector 

tf ) =C«(M 1) )C«(^ 1) )---C (1) (^ ) )|0>. (A.44) 

The result is 

A(1) ( A ) = — T74T + E( coth ( A + ^)+ coth (A-/i 4 )) 
smh(A) fr[ 

m 

- £ (coth(A + + coth(A - , (A.45) 
z=i 

where /4 satisfy the constraints 

n 

/W = ^(coth^f+^+coth^f-^)) 
1=1 

-2 £ (coth^+^+ooth^-^)) 
= 0. (A.46) 
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